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Abstract. In this work we present a flexible tool for tumor progression, 
which simulates the evolutionary dynamics of cancer. Tumor progression 
implements a multi-type branching process where the key parameters 
are the fitness landscape, the mutation rate, and the average time of 
cell division. The fitness of a cancer cell depends on the mutations it 
has accumulated. The input to our tool could be any fitness landscape, 
mutation rate, and cell division time, and the tool produces the growth 
dynamics and all relevant statistics. 

1 Introduction 

Cancer is a genetic disease which is driven by the somatic evolution of cells [1] , 
where driver mutations for cancer increase the reproductive rate of cells through 
different mechanisms, e.g. evading growth suppressors, sustaining proliferative 
signaling, or resisting cell death [2]. Tumors are initiated by some genetic event 
which increases the reproductive rate of previously normal cells. The evolution 
of cancer (malignant tumor) is a multi-step process where cells need to receive 
several mutations subsequently |3]. This phase of tumor progression is charac- 
terized by the uncontrolled growth of cells |2l4j . The requirement to accumulate 
multiple mutations over time explains the increased risk of cancer with age. 

There are several mathematical models to explain tumor progression and the 
age incidence of cancer |5I6I7| . The models have also provided quantitative in- 
sights in the evolution of resistance to cancer therapy [8] . The models for tumor 
progression are multi-type branching processes which represent an exponentially 
growing heterogeneous population of cells, where the key parameters for the pro- 
cess are: (i) the fitness landscape of the cells (which determine the reproductive 
rate), (ii) the mutation rate (which determines the accumulation of driver mu- 
tations), and (iii) the average cell division time (or the generation time for new 
cells). The fitness landscapes allow the analysis of the effects of interdependent 
(driver) mutations on the evolution of cancer [5]. 

In this work, we present a very flexible tool (namely, TTP, tool for tumor 
progression) to study the dynamics of tumor progression. The input to our tool 
could be any fitness landscape, mutation rate, and cell division time, and the tool 
generates the growth dynamics and all relevant statistics (such as the expected 



tumor detection time, or the expected appearance time of surviving mutants, 
etc). Our stochastic computer simulation is an efficient simulation of a multi- 
type branching process under all possible fitness landscapes, driver mutation 
rates, and cell division times. 

Our tool provides a quantitative framework to study the dynamics of tumor 
progression in different stages of tumor growth. Currently, the data to under- 
stand the effects of complex fitness landscapes can only be obtained from patients 
or animals suffering the disease. With our tool, playing with the parameters, 
once the real-world data is reproduced, the computer simulations can provide 
many simulation examples that would aid to understand these complex effects. 
Moreover, once the correct mathematical models for specific types of cancer are 
identified where the simulations match the real-world data, verification tools for 
probabilistic systems can be used to further analyze and understand the tumor 
progression process (such an approach has been followed in [TU] for the verifica- 
tion of biological models) . In this direction, results of specific fitness landscapes 
of our tool have already been used in a biological application paper [TT] . While 
we present our tool for the discrete-time process (which provides a good approx- 
imation of the continuous-time process) , results of our tool for the special case of 
a uniform fitness landscape in the continuous-time process have also been shown 
to have excellent agreement with the real-world data for the time to treatment 
failure for colorectal cancer |S]. 

2 Model 

Tumor progression is modeled as a discrete-time branching process (Galton- 
Watson process [32] )■ At each time step, a cell can either divide or die. The 
phenotype i of a cancerous cell determines its division probability bi and is 
encoded as a bit string of length four (i.e. {0, l} 4 ). The death probability di 
follows from 6j as d$ = 1 — bi. If a cell divides, one of the two daughter cells can 
receive an additional mutation (i.e., a bit flips from wildtype to the mutated 
type 1) with probability u in one of the wildtype positions (e.g., cells of phenotype 
1010 can receive an additional mutation only at positions two and four; cells 
of phenotype fill can not receive any additional mutations). The branching 
process is initiated by a single cell of phenotype i — 0000 (resident cell). The 
resident cells are wildtype at all four positions and have a strictly positive growth 
rate (i.e., 6 00 oo - doooo > 0). 

Fitness landscapes. Our tool provides two predefined fitness landscapes for 
driver mutations in tumor progression: (I) Multiplicative Fitness Landscape 
(Mfl) and (2) Path Fitness Landscape (Pfl). Additionally, the user can also 
define its own general fitness landscape (Gfl). A fitness landscape defines the 
birth probability bi for all possible phenotypes i. Following the convention of the 
standard modeling approaches, we let 60000 =1 / 2 (l + s o) be the birth probability 
of the resident cells (i.e., cells of phenotype 0000) |9lllj . The growth coefficient 
Sj indicates the selective advantage provided by an additional mutation at po- 
sition j in the phenotype. 



Multiplicative fitness landscape. In the Mfl a mutation at position j of the 
phenotype i of a cell results in a multiplication of its birth probability by (1 + s,-). 
Specifically, the birth probability bi of a cell with phenotype i is given by 

1 4 
&< = 2( 1 + s o)n( 1+? i)5 

3=1 

where = if the j-th position of i is 0; otherwise = Sj. Hence, each 
additional mutation can be weighted differently and provides a predefined effect 
(l + sx), (I + S2), (I + S3), or (I + S4) on the birth probability of a cell. Additional 
mutations can also be costly (or neutral) which can be modeled by a negative Sj 
(or sj = 0). If so — s i = s 2 = s 3 — s 4 the fitness landscape reduces to the model 
studied by Bozic et al. [5], which we call Emfl (Equal Multiplicative Fitness 
Landscape) and is also predefined in our tool. 

Path fitness landscape. The Pfl defines a certain path on which additional mu- 
tations need to occur to increase the birth probability of a cell. The predefined 
path can be 0000 -> 1000 -> 1100 -> 1110 -)■ 1111, and again the growth coef- 
ficients Sj determine the multiplicative effect of the new mutation on the birth 
probability (see appendix for more details). Mutations not on this path are dele- 
terious for the growth rate of a cell and its birth probability is set to 1/2(1 — v). 
The parameter v (0 < v < 1) specifies the disadvantage for cells of all phenotypes 
which do not belong to the given path. 

General fitness landscapes. Our tool allows to input any fitness landscape as 
follows: for bi for i £ {0, l} 4 , our tool can take as input the value of bi. In this 
way, any fitness landscape can be a parameter to the tool. 

Density limitation. In some situations, a tumor needs to overcome current 
geometric or metabolic constraints (e.g. when the tumor needs to develop blood 
vessels to provide enough oxygen and nutrients for further growth |13lllj ) . Such 
growth limitations are modeled by a density limit (carrying capacity) for vari- 
ous phenotypes. Hence, the cells of a phenotype i grow first exponentially but 
eventually reach a steady state around a given carrying capacity Ki. Only cells 
with another phenotype (additional mutation) can overcome the density limit. 
Logistic growth is modeled with variable growth coefficients Sj = Sj(l — Xi/Ki) 
where Xi is the current number of cells of phenotype % in the tumor. In this 
model, initially Sj ~ Sj (Xi -C Ki), however, if Xi is on the order of Ki, Sj 
becomes approximately zero (details are given in the appendix). 

3 Tool Implementation & Experimental Results 

Our tool provides an efficient implementation of a very general tumor progression 
model. Essentially, the tool implements the above defined branching processes 
to simulate the dynamics of tumor growth and to obtain statistics about the 
expected tumor detection time and the appearance of additional driver mutations 
during different stages of disease progression. TTP can be downloaded from 
here: http : / /pub . ist . ac . at/ttp[ 



For an efficient processing of the discrete-time branching process, the stochas- 
tic simulation samples from a multinomial distribution for each phenotype at 
each time step |9lllj . The sample returns the number of cells which divided 
with and without mutation and the number of cells which died in the current 
generation (see the appendix for details). From the samples for each phenotype 
the program calculates the phenotype distribution in the next generation. Hence, 
the program needs to store only the number of cells of each phenotype during 
the simulation. This efficient implementation of the branching process allows 
the tool to simulate many patients within a second and to obtain very good 
statistical results in a reasonable time frame. 
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Fig. 1. Experimental results illustrating the variability of tumor progression. In pan- 
els a and b we show examples for two particular simulation runs where the cells grow 
according to the Emfl and resident cells (blue) are constrained by a carrying capacity 
of i^oooo = 10 6 . In panel c the cells grow according to the Pfl. In panel d we show 
statistical results for the probability density of tumor detection when cells grow accord- 
ing to different fitness landscapes. Parameter values: growth coefficients so = 0.004, 
si = 0.006, s 2 = 0.008 , s 3 = 0.01, s 4 = 0.012, and v = 0.01, mutation rate u = 10" 6 , 
cell division time T = 3 days, tumor detection size 10 9 cells. 



Modes. The tool can run in the following two modes: individual or statistics. 
In the individual mode the tool produces the growth dynamics of one tumor in a 
patient (see panels a, b, and c in Fig. [I]). Furthermore, both the growth dynamics 
and the phcnotypc distribution of the tumor are depicted graphically. In the 
statistics mode the tool produces the probability distribution for the detection 
time of the tumor (see panel d in Fig. [T]) both graphically and quantitatively. 
Additionally, the tool calculates for all phenotypes the appearance times of the 
first surviving lineage, the existence probability, and the average number of cells 
at detection time. 

Features. TTP provides an intuitive graphical user interface to enter the pa- 
rameters of the model and shows plots of the dynamics during tumor progression, 
the phcnotype distribution, or the probability density of tumor detection. These 
plots can also be saved as files in various image formats. Furthermore, the tool 
can create data files (tab-separated values) of the tumor growth history and 
the probability distribution of tumor detection for any set of input parameters 
(details on the format are given in the appendix). 

Input parameters. In both modes, the tool takes the following input param- 
eters: (i) growth coefficients sq, si, S2, S3, and S4 (and v in the case of Pfl), 
(ii) mutation rate u, (iii) cell generation time T, (iv) fitness landscape (Mfl, 
Pfl, Emfl, or Gfl with the birth probability for each phenotype), and optional 
(v) density limits for some phenotypes. In the individual mode, additionally, the 
user needs to provide the number of generations which have to be simulated. In 
the statistics mode, the additional parameters are: the tumor detection size and 
the number of patients (tumors which survive the initial stochastic fluctuations) 
which have to be simulated. 

Experimental results. In panels a, b, and c of Fig. [I] we show examples of the 
growth dynamics during tumor progression. Although we used exactly the same 
parameters in panels a and b, we observe that the time from tumor initiation 
until detection can be very different. In panel d we show the probability density 
of tumor detection under various fitness landscapes. Further experimental results 
are given in the appendix. 

Case studies. Several results of these models have shown excellent agreement 
with different aspects of real- world data. In |9j, results for the expected tumor 
size at detection time using a Emfl fit the reported polyp sizes of the patients 
very well. Similarly, using a continuous-time branching process and a uniform 
fitness landscape, results for the expected time to the relapse of a tumor after 
start of treatment agree thoroughly with the observed times in patients [S]. 
Future work. In some ongoing work, we also investigate mathematical models 
for tumor dynamics occurring during cancer treatment modeled by a continuous- 
time branching process. Thus an interesting extension of our tool would be to 
model treatment as well. Another interesting direction is to model the seeding of 
metastasis during tumor progression and hence simulate a "full" patient rather 
than the primary tumor alone. Once faithful models of the evolution of cancer 
have been identified, verification tools such as PRISM JT3] and theoretical results 
such as [H] might contribute to the understanding of these processes. 
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A Appendix: Details of the Tool 



TTP is available for download at |http : //pub . ist . ac . at/ttp[ The tool is im- 
plemented in Java and runs on all operating systems which run a Java Virtual 
Machine (JVM) of version 1.7 or above ( |www . j ava . c om/get j ava ) . All the nec- 
essary libraries are included in the jar-file. 



A.l Tool Features 

Our tool supports various features in two running modes. 

In the individual mode, TTP simulates the tumor growth dynamics for a 
given number of generations. Plots of the growth dynamics over time and the 
current phenotype distribution are produced simultaneously. (Both plots can be 
saved in a PNG-file or SVG-file.) The full growth history for all cell types can 



also be stored in a data-file (format is described in section A. 2 1. 

In the statistics mode, TTP simulates the given number of patients with 
the same parameters and simultaneously shows the probability density of tumor 
detection for a given detection size (10 9 cells correspond to a tumor volume of 
approximately 1cm 3 ). The average tumor detection time and the average fraction 
of resident cells at detection are also shown during the simulations. After all 
patients have been simulated, the existence probability at detection, the average 
number of cells, and the average appearance year of the first surviving cell for 
all phenotypes are calculated and shown in a new window. In addition the tool 
shows the number of detected and died tumors per year in a separate window. 



All these data is stored in a data-file (format is described in section A. 2) 



A. 2 Installation and Implementation Details 

TTP is written in Java and makes use of several other libraries. The tool requires 
the Java Runtime Environment (JRE) of version 1.7 or above. To start TTP 
double-click on ttp.jar or on the command line type: java -jar ttp.jar. 
(Make sure that you have the permission to execute ttp.jar. On Mac OS in- 
voking the tool from the command line can overcome the security restrictions.) 

The tool is composed of the following components: the model implementation, 
the statistics thread, the graphical user interface, and the plot generator. 
Model implementation. The core component of the tool is the efficient im- 
plementation of the discrete-time branching process. Following Bozic et al. [9], 
the number of cells Xi of phenotype i in the next generation it + 1) is calculated 
by sampling from the multinomial distribution 

Probp^^) - (Vi,V»,V3)] = ^f^M 1 - «) Wl ^ (M)**] (i) 

where yi + y 2 + 2/3 = Xi(t), y 3 = J2k and 



Xi{t + 1) = Xi(t) + Y 1 - Y 2 + M kl 

k 



(2) 



The number of cells which give birth to an identical daughter cell is denoted 
by Yx, the number of cells which die is denoted by Y 2 . The number of cells 
which divide with an additional mutation is given by Y3 and the number of 
cells mutated from phenotype k to i is given by M^i- In general, one can define 
a mutation matrix to encode the probabilities rriki that a cell of phenotype 
k mutates to a cell of phenotype i. In our case this matrix is defined by the 
sequential accumulation of mutations. A cell of some phenotype can receive an 
additional mutation only on the positions in its bit-string encoding which are 
wildtype (i.e. only bit flips from to 1 are allowed). Mutations on all allowed 
positions are equally likely. Back mutations are not considered. 

Fitness landscapes. Our tool supports four fitness landscapes for additional driver 
and passenger mutations: (i) Mfl, (ii) Emfl, (iii) Pfl, and (iv) Gfl. In princi- 
pal, driver mutations increase the birth rate of a cell whereas passenger mutations 



have no effect on the cell's birth rate In the tables Al and A2 we present 
the complete definition of the Mfl and the Pfl, respectively. (The definition of 
Emfl and Gfl have been given in Section [2]) 



Table Al. Multiplicative Fitness Landscape. 



Additional mutations 


Phenotype 


Birth probability 





0000 


V2(l + so) 


1 


1000 
0100 
0010 
0001 


1/2(1 + S ) (1 + Sl) 

1/2(1 + s ) (l + s 2 ) 
1/2(1 +s )(l + S3) 
1/2(1+ So) (l + s 4 ) 


2 


1100 
1010 
1001 
0110 
0101 
0011 


1/2(1 + S )(1 + Si) (l + s 2 ) 
1/2(1 + S )(l + Sl)(l + S3) 
1/2(1 + So) (1 + Sl)(l + S 4 ) 
1/2(1 + So)(l + S 2 )(l + S 3 ) 
1/2(1 + S )(l + S 2 )(l + Si) 

1/2(1 + so)(l + s 3 )(l + s 4 ) 


3 


1110 
1101 
1011 
0111 


1/2(1 + S )(l + Sl)(l + S 2 )(l + S3) 
1/2(1 + S )(1 + S 1 )(1 + S2 )(1 + S4 ) 
1/2(1 + S )(l + Sl)(l + S 3 )(l + S 4 ) 
1/2(1 +S )(l+Sl)(l + S 3 )(l + S 4 ) 


4 


1111 


1/2(1 + s)(l + Sl)(l + s 2 )(l + s 3 )(l + Si) 



Density limit. Our tool allows a separate carrying capacity Ki for each pheno- 
type i. When the Gfl is used, in the beginning of the simulation, the growth 
coefficients Si are calculated from the given bi for all phenotypes i € {0, l} 4 since 
the density limiting effects are based on the values for Si. As a technical detail, 
for sizes Xi for which the birth probability bi would fall below (or for which, 
equivalently, Si would fall below -1), we set bi = 0. 

Statistics thread. The statistics thread handles the simulation of many iden- 
tical branching processes to obtain the statistical results. These simulations run 



Table A2. Path Fitness Landscape. 



Additional mutations 


Phenotype 


Birth probability 





0000 


1/2(1 + s) 




1000 


1 /2(1 + S )(1 + S 1 ) 




0100 


1/2(1 - v) 


1 


0010 


1/2(1 - v) 




0001 


1/2(1 - v) 




1100 


1/2(1 +s )(l + si)(l + s 2 ) 




1010 


1/2(1 - v) 


2 


1001 


1/2(1 - v) 




0110 


1/2(1 - «) 




0101 


1/2(1 - v) 




0011 


1/2(1 - u) 




1110 


1/2(1 + S )(l + Sl)(l + S2)(l + S 3 ) 




1101 


1/2(1 - v) 


3 


1011 


1/2(1 - v) 




0111 


1/2(1 - «) 


4 


1111 


1/2(1 + s )(l + Sl)(l + s 2 )(l + s 3 )(l + S A ) 



in a separate thread such that the GUI keeps responsive for user requests. After 
completing all the necessary simulations a data-file with all relevant results is 
automatically generated and stored to the execution directory of the tool. 

Graphical user interface. The graphical user interface (GUI) component con- 
tains frames and forms required for the functionality of the tool. It also handles 
all the user requests and distributes them to the other components. Within the 
GUI the plots for the tumor progression dynamics (in individual mode) and for 
the probability density of tumor detection (in statistics mode) are displayed. 
Multiple screenshots of the GUI are shown in section |A.4| 

Plot generator. The plot generation is based on the free JFreeChart library 
(www.jfree.orgl. For the generation of the Scalable Vector Graphics (SVG) the 
Apache XML Graphics library (http://xmlgraphics.apache.org) is used. An 
example for a plot generated by our tool is shown in Figure |A2| 

Data files. TTP produces various data files which can be used for further 
analysis and processing. The data are given as tab-separated values where each 
record is one line of the text files. In Listing [AT] we show an example of a data file 
generated in the statistics mode. Average results are given as comments which 
start with a hash. In the individual mode the data file contains the number of 
cells Xi for each phenotype i in all generations. 



Multiplicative Fitness Landscape - Total number of cells: 2.96E9 




0.0 2.5 5.0 7.5 10.0 12.5 15.0 17.5 20.0 22.5 25.0 27.5 

Tumor time (years) 



- 0000 


= 0001 


° 1000 


<- 0010 


0100 


a 1010 


^ 1001 1100 


" 0011 


° 0110 


° 0101 


• 1011 


' 1110 


« 1101 


* 0111 


□ 1111 





Fig. A2. Example for a generated plot of the tumor growth dynamics. 



Listing Al. Generated data file in the statistics mode. 

# Used fitness lands cap e : Multiplicative 

# Growth coefficients: s0 = 0.004, sl=0.006, s2 = 0.008, 88=0.01, s4=0.012 

# Mutation rate : IE— 6 

# G eneration time ( days ) : 3.0 



# Simulation f 


or 10000 


p atients . 








# Bin size 


i s 


20 generations. 










generation 


detected 


cumul 


_det 


died 


cumuLdied 















0.046003 





920069 


20 












0.001938 





958826 


40 












0.000705 





972923 


60 












0.000358 





98009 


80 












0.000227 





984622 


100 












0.000149 





987594 


5980 







1 







1 




6000 







1 







1 




# 0.7981% 


( abs 


: 10000) 


have been detected within 




at most 


# 


20000 generations 










# 99.2% (abs : 


1243006) 


tumors 


went 


extinct . 






# Average 


year 


of detection: 


22.04 









# Mutant appearance times : 

# Mutant 0000 appeared in average in g eneration : (year : ) , 

# exi stance probability: 100% , 

# average number of cells : 3267013.74 

# Mutant 0001 appeared in average in generation : 1540 ( year : 12.66), 

# existance pr ob ability : 100% , 

# average number of cells : 1925660.02 



# 1253006 runs have been p erform ed to create 10000 p atients 



A. 3 User Manual 



TTP is invoked by a double-click on ttp.jar or by the command java -jar 
ttp. jar. After the tool has started, the GUI can be used for all operations (see 



Figure A3 for a screenshot of the GUI) . 

Input parameters. In the control panel the tool takes all the main parameters 
for tumor progression: the fitness landscape, the mutation rate, and the cell 
division time. If one of the prespecified fitness landscapes (Mfl, Pfl, or Emfl) 
is used, the relevant growth coefficients have to be defined. If the general fitness 
landscape is used, a pop-up window appears after the selection of the Gfl and 
the specific birth probabilities for all the phenotypes can be defined. To add 
density limits for specific phenotypes, a pop-up window appears after "density 
limit" has been checked. For each phenotype a different density limit can be 
given (—1 indicates that there is no limitation on this phenotype). To obtain 
statistical results, the number of patients (i.e. the number of tumors with a 
surviving lineage) and the tumor detection size (i.e. the number of cells when a 
tumor can be detected) need to be provided. 

Modes. After all the parameter values have been specified, the tool can either 
run in the individual or the statistics mode. To simulate the growth dynamics 
of a single tumor, click on "New simulation" and the tool runs in the individual 
mode. Then any number of cell generations can be simulated until the tumor 
consists of more than 10 15 cells. 

The statistics mode can be started by clicking on "Obtain statistics". The 
tool simulates the given number of tumors until they reach the detection size 
and calculates all relevant statistics. 

Output. In the individual mode, the tool generates plots for the growth dy- 
namics and the phenotype distribution during the simulation. Furthermore, the 
entire tumor growth dynamics for each phenotype can be stored as a data file 
and the plots can be saved as PNG and SVG files (plots are stored to the folder 
"charts" in the execution directory of the tool). 

In the statistics mode, the tool generates the plot for the probability density 
of tumor detection. Statistics about the appearance time of the mutants and the 



detection and extinction year are shown in separate windows (see Figures A4 



and A6 for screenshots) . All the generated statistics are automatically saved to 



a data file (see Listing Al ) 



A. 4 Experimental Results Si Screenshots 

In this section we present some additional experimental results and multiple 
screenshots of the tool. In Table |A"3| we compare the probability of tumor detec- 
tion for the fitness landscapes Emfl, Mfl, and Pfl. In average our tool needs 
approximately 90ms to simulate a tumor with 10 9 cells (on a dual core 2.67GHz 
processor) . 



Table A3. Cumulative probability of tumor detection for different fitness landscapes. 
Results are averages over 10 runs. Parameter values: growth coefficients so = 0.004, 
si = 0.006, s 2 = 0.008, s 3 = 0.01, s 4 = 0.012, v = 0.006, mutation rate u = 10~ 6 , 
detection size M = 10 9 . 



generation 


Emfl 


Mfl 


Pfl 


2000 


0.0 


0.001 


0.0 


2200 


0.0 


0.038 


0.0 


2400 


0.0 


0.212 


0.006 


2600 


0.004 


0.491 


0.075 


2800 


0.068 


0.724 


0.288 


3000 


0.289 


0.865 


0.559 


3200 


0.567 


0.937 


0.766 


3400 


0.774 


0.971 


0.887 


3600 


0.892 


0.987 


0.947 


3800 


0.949 


0.994 


0.976 


4000 


0.977 


0.998 


0.989 


4200 


0.989 


0.999 


0.995 


4400 


0.995 


1.0 


0.998 


4600 


0.998 


1.0 


0.999 


4800 


0.999 


1.0 


1.0 


5000 


1.0 


1.0 


1.0 



1 1)1 Tool for Tumor Progression 



a 



Control 

growth coefficients 



0.004|jj s o | 0.006|j-| s i | 0.008[rj s 2 | 0.0l|-r| s 3 | 0.012|jj s 4 | 0.006|jj v 

fitness landscape 



other parameters 






mutation rate: 


0.000001 




generation time (days]: 


3 




density limit 







® path 

O multiplicative 
O general 



New simulation 



10p|-r| generations 



Simulate 100 generations 



Save PNG 



Save SVG 



Save Dynamics 



patients: 
detection size: 
generation: 3000 



1,000 








Obtain statistics 


1,000,000,000 









tumor size: 



4.99E12 



1 E 13 
1E12 

1 El 1 

1E10 

1EH 

1ES 

2 1E7 
m 

O 1EB 
IBS 
1E4 
1E3 
1E2 
1E1 
1E0 




0.0 2.5 5.0 7.5 10.0 12.5 15.0 17.5 20.0 22.5 

Tumor time (years) 



25.ii 



- 0000 = 0001 o 0010 
° 1011 * 1111 0011 



0100 t 1000 « 1010 t- 1100 □ 1001 - 1110 
0110 * 0101 □ 0111 



1101 



Fig. A3. Graphical User Interface of TTP in the individual mode. 
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Fig. A4. Graphical User Interface of TTP in the statistics mode. 
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Fig. A5. Statistical results of the average detection and extinction year. 
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Fig. A6. Statistical results of the average appearance year, the existence probability 
and the number of cells at detection time. 



